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Abstract 

Supersymmetric GUTs are the most natural extension of the Standard model 
unifying electroweak and strong forces. Despite their indubitable virtues, among 
these the gauge coupling unification and the quantization of the electric charge, one 
of their shortcomings is the large number of parameters used to describe the high 
energy thresholds (HET), which are hard to handle. We present a new method 
according to which the effects of the HET, in any GUT model, can be described by 
fewer parameters that are randomly produced from the original set of the parameters 
of the model. In this way, regions favoured by the experimental data are easier to 
locate, avoiding a detailed and time consuming exploration of the parameter space, 
which is multidimensional even in the most economic unifying schemes. To check 
the efficiency of this method, we directly apply it to a SUSY SO(10) GUT model 
in which the doublet-triplet splitting is realized through the Dimopoulos-Wilczek 
mechanism. We show that the demand of gauge coupling unification, in conjunction 
with precision data, locates regions of the parameter space in which values of the 
strong coupling a s t r0 ng are within the experimental limits, along with a suppressed 
nucleon decay, mediated by a higgsino driven dimension five operators, yielding 
lifetimes that are comfortably above the current experimental bounds. These regions 
open up for values of the SUSY breaking parameters mo,M 1 / 2 < 1 TeV being 
therefore accessible to LHC. 
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1 Introduction 



Grand unified theories (GUTs) provide a simple and elegant framework for the unification 
of strong, weak and electromagnetic forces. In addition, they offer a simple explanation of 
the electric charge and hypercharge assignments to the quarks and leptons in the Standard 
Model and combine its seekingly unrelated left and right - handed multiplets (five per 
family) into common representations of the larger unifying group. Moreover, their minimal 
supersymmetric versions, SUSY GUTs, lead to a successful gauge coupling unification 
at scales Mqut ~ 10 16 GeV, which is impossible to realize without supersymmetry. 
Also, in the framework of particulate SUSY-GUTs, the lightness of the neutrinos can be 
explained and a mechanism for Baryogenesis through thermal Leptogenesis is offered as 
an alternative to Baryogenesis through electroweak phase transition, which requires large 
CP - violating phases. 

In this note, our goal is to check the viability of SUSY GUTs, using electroweak 
precision and proton decay data by developing an integrated and simple scheme for the 
treatment of high energy thresholds (HET), which can be applicable, in principle, to any 
GUT model. In this scheme, the HETs are collectively parametrized by a small number 
of properly chosen variables, which are therefore easier to handle in phenomenological 
analyses. These are produced from the numerous parameters defining the HETs in a 
random way. 

We employ our method to a SUSY - SO (10) PQ, which seems to be a very promising 
candidate for a unified description. Models based on SO (10) unify all quarks and leptons 
of one family into one irreducible spinor representation, they give upper bounds on proton 
lifetime that still survive the current experimental bounds, they naturally incorporate 
the see-saw mechanism, reproducing the current neutrino oscillation data and explaining 
the lightness of the left-handed neutrinos with the sterile neutrino mass being of order 
Mqut 0, and finally predict Yukawa unification [3J, X t = \ b = \ T = A Pr , for large 
tan/3. Our analysis is based on a minimal supersymmetric 5*0(10) model first proposed 
in [1], whose low energy effective theory is the constrained MSSM. 

As we have already noted, a common characteristic of all SO (10) models is the in- 
clusion of all quarks and lepton of one generation, along with the right handed neutrino, 
in the same spinorial, 16, representation. The 5*0(10) generators belong to the adjoint 
representation 45y. 

As far as the content of the Higgs sector is concerned, there are two approaches 
that have been considered in the literature. Their difference rely on the way the Higgs 
mechanism is realized at the GUT scale. The first approach uses pairs of spinor Higgs 
multiplets in 16h + 16# representations to reduce the rank of the group [1H6], whereas 
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the other adopts a pair of 126// + 126// instead [7]. In both Higgs field in the 

45 or the 54 representation is also needed, in addition, to further break 50(10) to the 
Standard Model. Moreover, in order to build a viable model in 50(10) , one should always 
take into account that the Higgs multiplets of 50(10) should include the two MSSM 
Higgs doublets H u , Hd, which give masses to the up and down quarks respectively and 
additional Higgs fields, which are necessary to obtain the appropriate symmetry breaking 
pattern at the unification scale. The Higgs content of the particular model, which we will 
study, follows the first approach and is characterized as minimal. It consists of one Higgs 
multiplet A in the adjoint 45// representation, two pairs of 16// + 16// multiplets, named 
C + C ', C + C, and two Higgs T±, T 2 in the vector 10// representation. 

5O(10) is spontaneously broken to the Standard Model without intermediate break- 
ing scales, as shown in Figure [U 



5O(10) 




CTStandard ModeP> 



SU(3) C x SU{2) L x U(l) x U(l) 



B-L 



Figure 1: The breaking scheme of 50(10) to Standard Model. 



In general, there are at least four sectors that are needed in the Higgs superpotential 
to accomplish the 50(10) breaking [HE]: the doublet-triplet-splitting sector W 2 /3, the 
adjoint sector Wa, the spinor sector Wc and the adjoint-spinor coupling sector Wacc 
so that the total superpotential is 

W = W A + W c + Wacc + W 2/3 . 

The issue of the doublet-triplet splitting in 5O(10) is, in fact, a manifestation of the gauge 
hierarchy problem which is present in every grand unified theory. In SUSY 5O(10), 
the two Higgs doublets, H u , Hj, required for the electroweak symmetry breaking, are 
contained in the vector Higgs multiplets, H u C 5, Hd C 5 of SU(5), which fall into the 
10// representation of 5O(10). 

This pair of weak doublets must remain massless after the 5O(10) breaking at the 
GUT scale, while their color-triplet partners in the vector Higgs multiplet should obtain 
superheavy masses. The doublet-triplet splitting is implemented by the Dimopoulos - 
Wilczek mechanism [8]. This mechanism assumes the existence of a term of the form 
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T\ AT 2 in the Higgs superpotential and demands the adjoint Higgs multiplet A to have a 
vev along the B — L direction, 



The parameter a is of the order of the GUT breaking scale Mqvt- If one adds in the 
superpotential the mass term M 2 (T 2 ) 2 , half of the four Higgs doublets contained in T x , T 2 
remain massless and thus the picture of the MSSM Higgs spectrum is revealed at the low 
energy effective theory. Hence, the doublet-triplet-splitting sector will be assumed to have 
the following form, 



where A and M 2 are the massless and GUT scale massive parameters, of this sector, 
respectively. 

The adjoint sector, Wa, is responsible for the Dimopoulos - Wilczek vev (A) given 
in ([I]), which breaks SU(5) to the Standard Model symmetry preserving U(1)b-l ( see 
Figure Q] ) • 

The spinor sector, Wc, of the Higgs superpotential forces the pair of spinor Higgs 
multiplets C + C to get superheavy vevs along the SU(5>) singlet direction. In this way, 
the rank of the group is reduced from 5 to 4, since the B — L symmetry is broken. 

Finally, the adjoint-spinor coupling terms, Wacc'i are necessary to prevent the man- 
ifestation of colored pseudo-Goldstone bosons with small masses, which may destroy the 
unification of gauge coupling constants and the low energy particle spectrum. 

Apart from providing a decent Higgs mechanism, the model just described yields 
reliable results for the fermion masses [EJIH], including neutrino masses and their oscilla- 
tions [6l[9] j at the electroweak scale by introducing the proper Yukawa terms at the GUT 
scale. 

This paper is organized as follows. In Secj2j we briefly discuss the proton instability. 
The expected proton lifetime depends on a superheavy mass M e ff which can be expressed 
in an efficient way as a function of HETs and the gauge coupling constants at Mz- In 
Sectj3j we propose a new method to count in the effects of the HETs in the running 
of the gauge couplings, needed for finding M e ff , imposing the experimental constraints 
from a st rong and proton decay. The results of our analysis are presented in SecJH and our 
conclusions are given in SecJHJ 

2 Proton Decay 

In every GUT model, the baryon number violation is unavoidable and consequently nu- 
cleons decay. In 50(10) proton decays through D=5 operators, which are induced via 



(A) = diag(a, a, a, 0, 0) Cg) ir 2 . 



(1) 



W 2/3 = XTi A T 2 + M 2 Tl 



(2) 
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the exchange of superheavy color triplet Higgsinos in the 10 Higgs multiplet [10Tu3] . with 
the dominant decay mode being p — > V K + . These operators arise from an effective 
superpotential which is inverse proportional to a mass parameter to be denoted by 

M,„ = *M = <W (3) 
e// M 2 M 2 ' K J 

where M 3 M' z are the superheavy color triplet Higgs masses of the vector Higgs multiplets 
Ti, T 2 . Thus, the expected proton lifetime turns out to be proportional to M^. 
For the dominant decay mode, the decay rate is: 

T(p -> vK + ) = T (P -> "i K+ )- ( 4 ) 

i = e, fx, t 

Each of the partial widths in ([1J) are given by [T2HH]: 

Men, 

The factor f3 p denotes the hadronic matrix element between the proton and the vacuum 
state of the 3 quark operator [15], employed in the appropriate chiral Lagrangian schemes 
[To] . In our approach, we use [T7] 

P = (5.6 ±0.5) x 10" 3 GeV 3 

extracted from lattice gauge calculations. A in ([5]) depends on quark masses (at 1 GeV) 
[TTj and CKM matrix elements and is given by 



T(p^ u t K + ) =1^1 



A = m s m c Vl V 21 A s A L . (6) 



°2 

The first parameter, As, represents the short-range renormalization effects between GUT 
and SUSY breaking scales [UJHSj, while the second one, Al, accounts for the long-range 
renormalization effects between SUSY scale and 1 GeV [TTJ. Their values are As — 0.93, 
for m t = 175, [TS] and A L ~ 0.32 (2-loop result), 0. 

In (j5J), the BiS are the functions that describe the dressing of the loop diagrams and 
are given by the formula 

Bi 1 ^(p 2B2i + P 3 -^ B3i ), (7) 
sin 2/3 m s V^ \ m c V 21 V 2 2 J 

where is the contribution of the jth generation particles in the loop with 

B 3% = F{ui, dj, W) + {dj eA . (8) 
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The functions F in (J8J) contain the corresponding loop integrals [M], with i — 1,2, 3 and 
j = 2, 3, while P2 and P3 are inter-generational, CP violating phases given by 



e* 7 \ 



5>i = 0, i = 1,2,3. (9) 



There exist two distinct limiting cases having to do with the relative contributions of 
the second and third generation: the destructive interference, occurring for P3/P2 = — 1, 
and the constructive one when P3/P2 = +1. We adopt the second case to achieve 
maximum mixing, and hence smaller, lifetimes that are more tightly constrained by data. 

Finally, the factor C in (jSJ) contains chiral Lagrangian factors, which convert a La- 
grangian involving quark fields to the effective Lagrangian involving mesons and baryons 
|19j . Its value has been calculated to be C = 1.014 , according to the values given in [T4"] . 

The current experimental lower bound on proton lifetime from Super-Kamiokande 
is [20] 

T ip ^uK+) > 1-6 x 10 33 yrs. (10) 

For given SUSY inputs, this constrains the value of the parameter M e //, which, as we shall 
see, depends on the values of the gauge couplings at M z and other high energy threshold 
parameters of the theory. Note that in ref. [13], in the context of the SU(5) model, the 
color-triplet Higgs boson mass, which is the analogue of M e ff in 577(5), is constrained 
by the precision measurement bounds put on a stTong . Following an analogous treatment 



as in [13], we first solve the 1-loop RGEs for the gauge coupling constants CKj, in the DR 
scheme, from the GUT scale Mqut down to the electroweak scale M z and express M e ff in 
terms of aii(M z ). In this approach, we count in the high energy thresholds j2TJ[22] of the 
superheavy spectrum, as well as the low energy thresholds of all sparticles and heavy SM 
particles. The resulting expression, for the 50(10) is 

^ = /(*) , (11) 

M z 

with 

h(a~ l ) = ^ [3a^(M z ) - a^(M z ) - 2« 3 " 1 (M Z )] . 

The effect of the low energy thresholds in (11 II) are encoded within a» from the low energy 
boundary conditions that are imposed at M Z - In fact 

a x \M z ) = ^ a ^cos 2 ^(l-A 7 + ^ln Ai 



5 em wy 1 2tt M z 

a- 2 \M z ) = sin 2 Ml- A 7 + ^ In ^ ) 



a^(M z ) = a~l ong (M z ) |^ - -L + J_ l n ^ (12) 



z 
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where a em , a st rong are the electromagnetic and strong coupling constants respectively, 
8 W is the weak mixing angle and M$, Ms account for the low energy threshold corrections. 
For further details and for the definition of the remaining quantities in (Tl2|) see [23]. At 
a subsequent stage, one can include a correction factor to (11 ip to take care of the small 
two-loop corrections to the gauge coupling running. In our numerical analysis, the two- 
loop effects are properly counted for since we integrate the 2-loop RGEs of the gauge 
coupling constants in the DR scheme. More details on the method that we employ will 
be presented in the following section. 

The function f(x) in ffTTl) depends only on GUT physics details and, in particular, on 
the high energy thresholds of the superheavy particles involved. For the case of SO (10) , 
it is found that, 



/(*) 
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1 + 8x 2 \ 4 (1 + 4x 2 ) 3 



1/2 



(13) 



1QV2[\1+x 2 J 1 + 32x 2 
In this expression, the effect of the high energy thresholds of this particular model depend 
only on the massless, free parameter x, which is defined as 

* - h < 14 > 

with a, c being connected to the GUT scale vevs of the adjoint and spinor Higgs fields 
respectively. Obviously, the larger the x, the larger the function f(x) is, facilitating the 
satisfaction of the proton decay bounds. However, x is naturally expected to be of 0(1) 
as being the ratio of vevs which are both of order ~ Mqut • On these grounds, x cannot 
be taken arbitrarily large. 

For comparison, in the SU(5) model the function f(x) = 1 and we recover exactly 
the result presented in [13] . 



3 High Energy Thresholds 

It is clear from eq{TT]that, in order to derive M e ff , we need the values of gauge couplings 
at the electroweak scale. Their values are set by ffT2"]) in terms of Ms, Ms . These are not 
physical masses but rather a convenient device which encodes all information for the low 
energy supersymmetric thresholds and heavy standard model states. Given the SUSY 
inputs and a em , a s , sin 2 #„,, we solve the RGEs imposing the gauge coupling unification 
condition, 

ai(M GC7T ) = a 2 {M GUT ) = a 3 (M GUT ) = a G , (15) 

but we don't restrict further our analysis by insisting on Yukawa unification. At Mqut , we 
also impose universal boundary conditions induced by gravity mediated SUSY-breaking 
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for the soft supersymmetry breaking parameters, although other schemes may be available, 



rrii = m , M* = M 1/2 , A { = A . 

Using standard procedure, we calculate the gauge coupling constants at any scale \i 
below the SUSY thresholds by solving the appropriate RGEs incorporating the effects of 
the thresholds of all heavy Standard Model particles of mass msMi > A* as wen as those 
of all low energy SUSY particles Si and superheavy particles Hi with masses ~ Mqut 
associated with the specific GUT model at hand. The result is 

= Qg^^Mqut) + (2 — loops effects) 

+ ^ (bf M + bf USY ) In I*™*- + J- V b? M * In -JL. (16) 
2tt w 1 a 2 it ^ 1 m SM . 



GUT 



bf are the beta function coefficients of any species A with bf M , bf USY being the beta 
function coefficients of «j of the SM and low energy supersymmetric modes respectively. 

From the evolution of a* from Mqut to the lowest high energy threshold, Ml, one can 
easily derive, ignoring momentarily the two loop effects, that 



of ' W = a -a\M GVT ) + 1 6- ta^ff + ^E a? h. ^ , (17) 

Hi 

where 6 GC/T = b SM + b SUSY + b H . By defining 



1 - GUT i M GUT 



ocq 1 (M l ) = a G \M GUT ) + — b GUT In 



2vr Mi ' 

which represents the running from Mqut^o Ml, if HETs are ignored, and by using 

Q = — V 6f« ln^- (19) 

equation f lT7|) can be cast in the form, 

ar\M L ) = a G \M L ) + Q . (20) 

where the effect of HET is included within q. 

Eqs. ( 120]) can serve as boundary condition at the lowest HET, Ml that takes into 
account the effects of all HETs, which are included within the constant c, . Between Ml 
and the electroweak scale \x = Mz no high energy thresholds are present and RGEs run 
as usual. 
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The importance of the adoption of the parametrization ffT9]) lies on the fact that the 
three q's carry all information on the the masses of the non-singlet superheavy fields of the 
model which contribute to the HETs in the running from Mz to Ml and vice versa. For 
any set of the model parameters, say pj, we assign a "vector" in a five dimensional space 
c— (ci, c 2 , C3, Ml, Mqut), which includes, except q, the values of the maximum, Mqut , 
and lowest high energy mass, Ml ■ For instance, in the version of the 5*0(10) model we 
are considering, the number of the parameters pj is ten and for any point in this ten- 
dimensional parameter space there correspond twenty-five correlated superheavy masses 
from which we determine Mqut-, Ml and Cj s through their definitions in (fT§j) . Further, 
in order to utilize q as inputs, we use a random sample generator, which assigns random 
numbers to the GUT parameters pj. In this way, random points p = (pi,P2, ■■■Pn) are 
drawn in the model parameter space and each of this is mapped to a c defined before. 
Consequently, our analysis is fully constrained from the random sample results and instead 
of dealing with a large number of GUT parameters and masses, we only have a few to 
consider in our analysis, namely c\, C2, C3, Ml, Mqut, which define c. The random 
procedure actually makes a selection by mapping the parameter space to a rather confined 
region, at least in the SO(10) model, which is spanned by the vectors c . Then, within this 
region, points satisfying the experimental criteria can be sought and the region shrinks 
even more. Consequently, one avoids time-consuming scans over a multidimensional (10- 
dimensional in S'O(IO) ) parameter space, since the random procedure has already selected 
the points c which meet the criteria. In this way we have found a very convenient way 
to parametrize the effect of HET, which is applicable to almost any GUT model, using 
the variables c . 

We will now discuss in more detail the numerical procedure we follow. For any input 
point in the space of the randomly generated vectors, c in , we pick up the points cf\ c™ 
and by running the 2- loop RGEs upwards, starting from from Mz, we determine the 
values of M L and ac(M L ) where the boundary condition f l2"Uj) are satisfied. By f fTgj) the 
value of OiG^McuT ) is also determined. Note that in this way the input value, M™, is 
not the same with the extracted value, Ml ■ Subsequently, we define C3 so that 



and relation ( 12 Op is satisfied for the coupling a 3 . That done, the initial point c j n , with 
c- n , M[ n , has moved to another with coordinates c[ n , c™, C3, Ml, M^j T , i.e. the values 
of C3 and Ml have been only changed. This procedure is repeated in each iterative step 
from the electroweak to the GUT scale, in the usual manner, until convergence has been 
obtained. In each iteration step, c\, are also corrected to meet the unification criteria 
but these corrections are small. Thus, at the end, we get a final point, Cf in , which is a 




(21) 



S 



successful point if it belongs to the set of the randomly generated vectors c or unsuccessful 
and hence discarded if it lies outside the region spanned by c points. In order to test 
the correctness of our method, we check, at the end, if a successful point maps to itself, 
modulo small differences due to numerical accuracies. This procedure will be explained 
and quantified in more detail later. 

Obviously, the set of successful points that pass the unification test will be reduced if 
additional physics constraints are imposed, like, for instance, bounds put by proton decay. 

The advantage of the method is that the running, from Mz to Ml, of the couplings 
involved is done without the effect of the HET in the RGEs. Their effect is taken into 
account by the boundary condition (1201) . which incorporates all the HET information 
within Cj's. 

The boundary conditions for the couplings and soft mass parameters are imposed at 
the unification scale, and, in the running from Mqut to Ml, the HET play an essential 
role. Since, in our approach, we want to treat HET collectively, through quantities similar 
to Ci , without knowledge of their precise values, we use the two-loop RGEs without the 
contribution of the superheavy particles, which are known. However, in order to include 
the effect of the HETs, at the end we correct each derived quantity at the scale Ml as 

F cor (M L ) = F(M L ) + A F . 

The added quantity, A F , depends on q , on the beta function coefficients bf = ^2 Hi bf* 
of the superheavy modes that are known and on In (Mgut/Ml) ■ The correction Ap is 
different for each quantity F and, for its knowledge, the one-loop explicit dependence of 
the RGE of the F need to be known. In particular, if 

dF 



dlnQ 



^2 G i a i + ■- 



where the ellipses denote contributions not explicitly dependent on the gauge couplings, 
the correction is 

A _ ( ^G_ i„ ( M GU'. 

F 4tt 



This is a valid approximation provided that the lowest, Ml, and the highest, Mgut, 
threshold are not far apart. In particular, < 10~ 3 is demanded. From the random 

samples, we find that on an average log A ,^ T — 2.7 and thus the approximation is more 
than satisfactory. 

With these in mind, we calculate couplings and masses at the electroweak scale in the 
ordinary manner with the values of the gauge couplings determined by the electroweak pre- 
cision measurements, [2"U], namely the effective mixing angle. sirij9(Mz)(MS) = 0.23152(14), 
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the value of the strong coupling constant: a s (Mz) = 0.1176(20) and the electromagnetic 

2 A 

coupling a em . Note that the experimentally measured effective mixing angle siritd is not 
the same with sin 2 9 W appearing in (112]) . which is defined as ratio of couplings in the DR 
scheme. Actually, the two are related by 

sin ft = sin 2 9 W (1 + Afcj ) , 

where Akf is calculated by the effective Zff coupling. 

The minimization conditions are solved with all one-loop effective potential corrections 
and the dominant two-loop QCD and top Yukawa corrections taken into account. The 
value of the parameter is then determined by the minimization conditions and tan/3 
is input. Therefore, apart from the particular treatment of HETs, the procedure is the 
standard one encountered in the constrained MSSM models. As far as proton decay is 
concerned, for given SUSY inputs, m , M 1 / 2 , A , tan/3, we derive M e ff, from f lTT]) . and 
-Bj's, in eq.((7]), which both affect the proton decay. Then proton lifetime is calculated and 
it provides an additional constraint. Since the dependence of M e ff on HETs is explicit 
only on the quantity x , defined by eq. fTH|) . to facilitate the analyses we shall pick slices 
of fixed values of x , in the space of random points, within which M e // is almost constant. 

4 Numerical Analysis - Results 

We follow the procedure described in Secf2]and|nj in order to delineate the acceptable pa- 
rameter region of the model at hand, which complies with electroweak precision data and 
proton decay constraint. Satisfaction of the experimental bounds on the strong coupling 
constant a str0 ng , by itself, imposes severe constraints, as we shall see, in conjunction with 
precision measurements and unification conditions. 

In order to facilitate the analysis, we generate random samples of randomly generated 
points c for which Mqvt is fixed. In our analysis, we present results for Mqut = 
2 • 10 16 GeV but higher values are not excluded, yielding qualitatively similar results. 
However, perturbativity limits on Yukawa and gauge couplings poses upper bounds on 
higher M GUT values and hence such large values are not considered. Besides, we select 
points defining "slices" in the space of c vectors for which the ratio x defined in eq. (fl4l) is 
x = 5 . In fact, we have found that a larger x satisfies the proton decay constraint easier, 
while a smaller one fails on both proton decay and a s t r0 ng constraints. In the following, 
for a more clearer presentations, only 1000 points are displayed in each Figure. 

In general, the effective weak - mixing angle, denoted in the Figures by Sf, takes values 
with error less than 3cr over all the parameter space, but the strong coupling constant 
prefers rather low values of SUSY breaking parameters. This is expected since for high 
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values SUSY is absent, due to decoupling, and, in this case, we deal with a conventional 
GUT model for which gauge coupling constant unification is hard to achieve. On the other 
hand, high SUSY breaking parameters and, in particular Mx/2 , which affects wino masses 
and in turn Eq. (151), shrink the range of the allowed values of M e ff, by proton decay 
bounds, leaving out a small number of the initial randomly generated points that pass 
successfully the proton decay and precision data tests. This is the reason only relatively 
small values for the SUSY breaking parameters are considered, m , M\i 2 < 1.5 TeV . 

4.1 The Q s t,. OJls and proton decay constraints 

We come to the point of examining the dependence of the strong coupling constant 
& strong on the supersyniHietric and HET parameters. We first observe its variation with 
changing m and Mi i 2 ■ In Figure [210, we display the pairs of C\ , c 2 as they are randomly 
generated. The value of the unification scale has been taken Mqvt = 2 ■ 10 16 GeV and 
on naturalness reasons all high energy parameters having dimension of mass are randomly 
generated with values differing from Mqut by at most three orders of magnitude. Al- 
though these parameters have been randomly generated from the independent parameters 
of the model, in the way prescribed earlier, they are correlated as shown clearly in the 
Figures. Therefore, successful points ought to be within the diagonal stripe displayed in 
these Figures. 

The black points represent those pairs of c\, c 2 that fail even to give unification at the 
quoted Mqjjt , after the 2 - loop running of the RGEs. For the green points, unification 
has been achieved but the value of a stron g is more than 4<r away. The magenta points yield 
otstrong with error smaller than 4er, and the yellow region is the subset of magenta points 
corresponding to values of a str0 ng with the smallest possible error < 2a. 

In Figures |2] to [5J we vary the input values of mo and Mi/ 2 , as shown on each Figure, 
by increasing either one, or both, of m , Mi/ 2 , but keeping Mqut , tan/3 and A Q fixed. 
In Figure El we display the analog of Figure |2] but with c 2 replaced by C3. A correlation 
of Ci and C3 is also observed with points that are more widely scattered. 

By raising the values of mo and Mi/ 2 , we have observed a considerable decrease of the 
number of random sample points that result in an acceptable a. s trong value. The yellow 
region shrinks and shifts slightly towards smaller values of c 2 and eventually disappears, 
as is the case in the displayed Figures [3] to [5j In those cases, the deviation of the theo- 
retical from the experimental value of a strong never becomes less than 2a . In Figure HJ by 
increasing the value of mo, keeping Mi/ 2 = 800 GeV as in Figure El the yellow and the ma- 
genta regions disappear, which means that the predicted theoretical values for a strong are 
outside the experimental limits. Therefore, to be within experimental limits, the values 
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Figure 2: In this Figure tan /3 = 10 and A = 100 GeV, with m , M 1 / 2 and Mqut as displayed. 
For a detailed description of the Figure see main text. 




of M1/2 and primarily mo should be kept central to small. The a s t r0 ng constraint is best 
met in a region of the m - M1/2 plane bounded by values of mo up to 1000 GeV and Mi/ 2 
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Figure 6: As in Figure [21 



up to 1300 GeV. For larger values of mo up to 1400 GeV, we should choose Mi/ 2 between 
900 and 1300 GeV at most and for larger values of m up to 1500 GeV, we should go for 
small mo around 500 GeV. This conclusion was rather expected since for large values of 
m and Mu2, exceeding 1 TeV, the MSSM and the Standard Model, due to decoupling, 
give almost the same predictions for the unification of coupling constants, where SM fails 
to yield a satisfactory unified picture. 

Based on the most recent experimental bound on proton lifetime (fTUj) . the lower bound 
on its lifetime translates into a lower bound on the mass parameter M e //that controls 
the proton decay width. This, denoted by M e ff(exp) , can be extracted from eq. §5§ and 
it is defined below. We perform these calculation for every randomly generated pair of 
ci, c 2 and the constraint of proton decay is satisfied provided that 

M eff (th) > M eff (exp) ee (3 P \A\ l 5 ^' ( 22 ) 

where is the bound in ([TO)) and M e ff(th) is read from (fTTl) . 

Running the RGEs we find points for which (]22j) is satisfied and the resulting M e fj(th) yields 
a proton lifetime in the range 10 34 — 10 37 years. In Figure [3 we illustrate this constraint 
for the case of mo = M1/2 = 800 GeV, by coloring in gray the pairs of c\, c 2 which fulfill 
(]22j) and in black the rest of those. These points overlap with the majority of the points 
that yield gauge coupling unification with values of a str0 ng within the la experimental 
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range. Raising of mo and/or M1/2 results to an homogeneous reduction of the gray region 
around a central point. Besides, the unilateral increase of mo keeping Mu% rather small 
leads to experimentally unacceptable proton decay rates. 




Figure 7: In this Figure tan/3 = 10 and A = 100 GeV, with m , Mi/ 2 , Mqut as quoted. For 
a detailed description of the Figure see main text. 

The results depend on the unification scale Mqut , which is an input in our analysis. 
For a str ong , raising M GUT by almost one order of magnitude ( M GUT = 9 • 10 16 GeV ) 
broadens the yellow region of Figure |2] which now includes smaller values for c±,C2- We 
note that almost half of the yellow points yield a strong within la error. Taking the same 
high value for Mqut , i n the case of the proton decay constraint, results to the same 
behavior. In fact, the gray region of Figure [7] moves towards smaller values of C2, while 
it reaches a wider range of values for c\. It is interesting to note that half of the points 
which bring unification of the coupling constants fulfill the proton decay constraint for 
Mgut = 9 • 10 16 GeV, while for M GUT = 2 • 10 16 GeV the corresponding rate is 40%. 
Thus, by pushing Mqut to higher value provides easier satisfaction of the constraints. 
As far as A is concerned, this parameter does not seem to play a significant role in our 
analysis. 

Our findings depend also on the value of tan/3. A change of tan/3 from 10 to 45 
causes a small decrease, ranging from 5 to 15%, in the number of points which succeed 
to give unification and, at the same time, shrinks their range of output values. A likely 
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explanation for this is that a high value of tan (3 increases the chance of a Landau pole to 
appear during the running of the RGEs. Also, the raise of tan /3 bounds the points which 
give a S f ron9 with error more than 4<r to have only positive c\ and only negative C3. On 
the other hand, with a large tan /3 the number of points which give a. strong with error less 
than 4a (or 2a for mo = M1/2 = 800 GeV) slightly increases. As far as the proton decay 
constraint is concerned, the points which satisfy f[2"2"j) show a considerable decrease, which 
starts from 22% for m = M1/2 = 800 GeV and reaches a 100% for m = 1500 GeV and 
M1/2 = 800 GeV. This was rather expected since Bi in (j2"2"l) depends on Hence, 
a change of tan/3 from 10 to 45 quintuples or so the values of M e ff(exp) leaving, at the 
same time, the values of M e ff(th) almost unchanged. Therefore, the chance of satisfying 
(j22p decreases. Overall, a change of tan (3 from 10 to 45 makes the satisfaction a harder 
task. 

4.2 Convergence of results 

In section [31 we pointed out that, in general, there are differences between the final points 
Cfi n = {c{ m , C2 , c{ m , M[ m }, stemming from the running of the RGEs (output values), 
and the initial (input values) c in = {c™, c 2 n , cf 1 , M|™}, generated from the random sample, 
owning to the corrections imposed to achieve unification at the given Mqut ■ For the case 
of Ci and c 2 , there is a close agreement between the input and output values but the same 
cannot be asserted for C3 as discussed in section [3j Actually, this deviation is expected 
since the C3 parameter is strongly correlated to a. strong and is forced to make the strong 
coupling be compatible with the unification scale that is determined by the couplings 0^2 • 
This is implemented by the shift of eq. (I21j) . resulting to c{ m , which is always towards 
higher values comparing to input c 3 n . This is also the case for C]" 2 and M™ , which are 
also shifted towards higher values but to a much lesser extend. 

Focusing only on the ability of the random sample points to accomplish unification 
of the gauge couplings, without any other experimental constraint, an overlapping region 
of both input and output points must exist for the model to be successful. We have 
found that the gap between the initial and final points augments by increasing mo and 
by reducing Mqut ■ 

The situation becomes more constrained if we simultaneously demand satisfaction of 
both Oistrong and proton decay constraints. For instance in Figure [8] we display in gray the 
points of C\j C3 as they are randomly generated. Actually, this represents a magnification 
of the central and upper area of Figure El The green points represent the output values 
c i ) c 3 °f tli e gray input points c^", c| n which yield values of a s4rong with error higher 
than 4cr. The magenta points are the corresponding points which yield a s t r0 ng with smaller 
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Figure 8: tan/3 = 10 and A = 100 GeV, with mo, M 1 / 2 , Mqut as quoted. For a detailed 
description of the Figure see main text. 



error, between 4a and 2cr, while the yellow points give a strong with error less than 2a. The 
gray crosses designate those output points that satisfy the proton decay constraint. 

We perform a "second" run using as inputs the elements of the previously produced 
output groups. In this run, we first check whether these points are mapped to themselves 
as they should. These are likely to be successful points provided they are subset of the 
initial randomly generated points. Obviously, only a subset of those survives, if it does at 
all, and this comprises the set of successful points. In order to check numerically whether 
a final point belongs to the set of the initial points, we define the "distance" Xi °f the 
point Cfi n from any c- n of the randomly generated points 



Xi = 



fin in, i 




fin in,i 
C 2 — C 2 


+ ... 


fin in, i 
C 5 - C 5 


in, i 


+ 


in, i 
C 2 


in, i 
C 5 



The minimum of these x = m ^ n {x'i s } defines the distance of c/j n from the set of the 
randomly generated points and if zero it means that the final point coincides with one of 
the random points that were initially created. In practice, the smallness of x indicates 
that the particular point is successful in the sense that it belongs to the 5*0(10) model 
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and in addition it agrees with the experimental results. 

In our analysis and with one million random points, we have found that points with 
X < 0.1 are acceptable by the model. 

This analysis certainly depends on the SUSY inputs. In the m - Mx/2 plane, successful 
points are found for values of mo and M\n reaching roughly 1000 GeV. If one keeps the 
higher end of M1/2 values constant, the restrictions are satisfied altogether for values of 
m up to approximately 1200 GeV. If we loosen the restriction for x to x < 0.2, keeping 
the other two constraints untouched, our successful region matches the one described in 
subsection 14.11 This findings supports not only our method but also the S'O(IO) model 
we have followed on the ground of satisfying a s t r0 ng and proton decay constraints. 

5 Conclusions 

We have presented a new approach towards treating in a collective way the large number 
of HET encountered in GUTs quantities by using a random sample technique. According 
to this technique the large number of the GUTs parameters, defining the high energy 
sector and hence the threshold masses, are mapped to a few properly defined parameters 
encoding all the information associated with the renormalization group running of the 
various quantities involved from the unification scale down to electroweak energies. This 
method is simple and efficient since : 

1. It avoids unnecessary runnings, occurring in the conventional scheme, where for 
each point in the multidimensional parameter space one has to solve numerically 
the RGEs. 

2. The parameter space is mapped to properly defined quantities, associated with 
the running of the gauge couplings, which collectively include enough information 
through which regions favoured by precision data on gauge couplings are easily 
located. 

3. The RG equations are solved without the inclusion of High Energy Thresholds 
and their effect is duly taken into account by changing appropriately the boundary 
conditions at the lowest of the high energy thresholds and at the Unification scale. 

4. This scheme is fast and accurate if the lowest, Ml, and the highest, Mqut, of 
the high energy threshold are separated by at most three orders of magnitude, 
M L /M GUT < 10- 3 . 
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We applied this technique to an 50(10) -based model where only five parameters 
( c i,2,3 tan# and Ml) embody all information associated with the HETs. Our analysis 
shows that this method is both convenient and efficient. We test it and demonstrate 
that there exists a confined region of values which yields results in agreement with 
the limitations imposed by precision data and coupling constant unification, as well as 
proton decay constraints. These regions are favored by central to small values of m and 
M1/2 in the the region of 500 GeV —1.5 TeV and are suppressed for large tan/3 due 
mainly to proton decay constraints. The effect of the common trilinear coupling Aq is 
small provided it lies in the TeV range. 
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